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SUMMARY 


The effects of solar radiation pressure on the 
motion of an artificial satellite are obtained, includ- 
ing the effects of the intermittent acceleration which 
results from the eclipsing of the satellite by the 
earth. Vectorial methods have been utilized to obtain 
the nonlinear equations describing the motion, and 
the method of Kryloff-Bogoliuboff has been applied in 
their solution. 
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THE EFFECT OF SOLAR RADIATION PRESSURE 


ON THE MOTION OF AN ARTIFICIAL SATELLITE 


by 


Robert W. Bryant 
Goddard Space Flight Center 


INTRODUCTION 


In the motion of artificial satellites with large area-to-mass ratios, such 
at the Echo I balloon satellite, the magnitude of the solar radiation pressure 
effect is substantial and must be taken into account in the analysis of the track- 
ing data. It is of fundamental importance in any study attempting to determine 
atmospheric density from satellite drag data to distinguish the energy loss (or 
gain) produced by radiation pressure from the loss due to drag. Furthermore, 
the effect is important since the orbital lifetime of a satellite will be substan- 
tially affected by the changes in energy resulting from radiation pressure. Sev- 
eral authors have considered the effects of solar radiation pressure on the motion 
of an artificial satellite (References 1,2, and 3); it is the purpose of this paper 
to present these effects, including that of the earth's shadow, in the form of non- 
linear equations derived vectorially and suitable for use with electronic computers. 


THE KRYLOFF-BOGOLIUBOFF METHOD 


H. Lass and J. Lorell (Reference 4) have applied the method of Kryloff- 
Bogoliuboff (Reference 5) in the following low-thrust cases: constant radial, 

constant transverse, and intermittent. 

The method may be described simply: A nonlinear equation 
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— = ^f(x, sin t) , 

where fi is a small number, becomes, upon integration. 


x ( t + 277 ) - x(t) 

277 


277 


f 


t + 2 t 7 


f(: 


sin r) dr . 


(1) 


Since y, is small, x may be considered constant during the integration from 
r - t to t — t + 2 tt. Since sin t is periodic, Equation 1 becomes 
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x(t+ 277)-x(t) y [ r/ \ J t-w \ / 

2n ^ 2 t 7 f(x, sin r) dr = /xF(x) . \£) 

Replacing the left side of Equation 2 with dx/dt, since the slope of the secant 
line is approximately that of the tangent, gives 
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j f ( x , sin t) dr = yF ( x ) . 


This result may be obtained simply by averaging the original equation over 
one cycle of the motion with respect to time. 


THE CHANGE IN ENERGY PRODUCED BY 
RADIATION PRESSURE 


The effect of radiation pressure, if the earth's shadow is neglected, 
will provide only short period terms in the development for the semi-major 
axis of the orbit (Reference 6). However, the inclusion of the shadow effect 
results in quite a different phenomenon. 

From the energy integral and Brown's 8/dt operator (Reference 6), we 
obtain 

2a 2 v • F 


da 

dt 


(3) 
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The velocity vector v may be represented in terms of a unit vector P 
directed from the center of the earth to perigee and a vector Q at right 
angles to P such that P, Q, R, with R the unit vector in the direction of 
the angular momentum vector, form a right handed system. The vector F 
is the disturbing acceleration, a is the semi-major axis, and /x = GM is the 
gravitational constant multiplied by the mass of the earth. 

We have 


v 


Qv//xa(l— e 2 ) COS E “ P \/fjL a sin E 
r 


and 


F = “ Fu° , 


where 


e = the eccentricity, 

E = the eccentric anomaly, 

F = the magnitude of the disturbing acceleration, 
u° - the unit vector in the direction of the Sun, and 


r = a ( 1 - e cos E) ; 
and we have 


P 


- cos i sin co sin fi + cos co cos fi 
+ cos i sin co cos fi + cos co sin fi 
+ sin i sin co 


0 


R 


— cos i cos co sin fi — sm co cos fi 
+ cos i cos co cos fl — sin co sin 0 
+ sin i cos co 


+ sin i sin fi 
- sin i cos fi 


+ cos i 
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Hence we obtain for Equation 3: 


da 

dt 


- 2a 2 F 
r V/xa" 


• u° a \/ 1 - e 2 


cos E 



(4) 


Considering the slowly varying quantities to be constant during one revolu- 
tion, and using 



The quantities Ej and E 2 are the eccentric anomaly at the exit and entrance 
of the shadow, respectively, and are obtained from the equation 

r • u° = - yjv 2 - p 2 ( 6 ) 


or, in terms of E, 


P • u® (cos E - e) + Q • u° yfl - e 2 sin E = “r ( 1 - e cos E) 2 - j , (7) 

where p is the radius of the earth. The solution of Equation 7 is readily 
obtained by the iterative technique and the entrance and exit roots are easily 
identifiable. We may write the expression for Aa as follows: 

Aa = - ^ sin — - f El V(Q • u°) 2 (l-e 2 ) + (P • u 0 ) 2 sin 
where 




Q • u° >/ 1 - e 2 

-P ” u° 


arctan 
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An examination of this expression shows that: (1) if there is no shadowing, 

sin (E 2 -E 1 )/2 = 0, and therefore Aa = 0; and (2) if e = 0, then <p becomes 
- (Ej + E 2 )/2 from a consideration of Equation 7, and again Aa = 0. 


THE COMPLETE SYSTEM OF EQUATIONS 
FOR THE OSCULATING ELEMENTS 


The equations for the elements may be obtained (Reference 7) as 
follows : 


where 


da 2a 2 v * F 


dco 

dt 


+ cos i 


d t 


de 

dt 


di 

rW 

dt 

yjfia ( 1 - e 2 ) 

dfi 

rW 

dt " 

\4ia( 1 - e 2 ) 

dfi 



cos (co + f) ; 


sin (co + f ) ; 




T = F R x r° 
N = F • Q , 

W = F • R , 


f = the true anomaly, 
L = F • P, and 


r° = the unit vector in the direction of r. 
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The complete procedure will be carried out for only one equation; the others 
may be obtained in exactly the same way: 

^ = yE (1 _ e 2 ) (T cos E + N) , (9) 

where 


T = 

- Fu° 

Rx r ® = " Fu ° RX r 

r 


- Fa 

Q * u° (cos E - e) - P * u® y/l - e 2 


r 

N = 

- Fu® 

• Q = -F Q * u® ( 1 - e cos E) , 

r 


r = Pa (cos E - e) + Qa \J\ - sin E . 


Substituting the above values into Equation 9, we have 


at = - VJ F [« • *° 7 (1 - 2 ' ” E ♦ I c. 2E) 

- P • u® “ \/l “ sin E cos E J ; (10) 


from which 


e( t + r) - e( t ) 


t + r p / 

■ -I J 


2e cos E + — cos 2E 

^ i 


- P • u° ^ v/]T 


e z sin E cos E 


dt , 


where r is the period of the satellite. 

Since F is small we may consider slowly varying quantities constant 
during the integration from t to t + r. By changing the independent vari- 
able for the integration and replacing the secant by the tangent on the left 
side of the above equation we obtain 
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if = f [”«■"• r*(f- 


2e cos E + ~2 cos 2E ] dE 


, r E 2 

- P • u® \ 1 - I sin E cos E dE 


(ID 


The limits are E^ to E 2 since F is zero outside this range. Therefore 
Equation 11 becomes 


de 

dt 


= “ ^ ( 1 - e 2 ) F 0 • u° E - 2e sin E + ^ sin 2E^ 

P ■ u° v/l - e 2 


The complete set of equations for the osculating elements is: 

"1 E . 
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do) 


dfi 


dF + cos i dt 


\Za~ F 
277 y/jT < 


«• (■ 


e cos E — ^ cos 2El 


- p • 


u° \/l - e 2 ^ 


3E . p 1 • 0If l 
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2 4 j 


J E, 


The integration of this system of equations can now be carried out with a 
large scale computer. 
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